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Abstract: We introduce a composition of quantum states of a bipartite system which is based 
on the reshuffling of density matrices. This non-Abelian product is associative and stems from 
the composition of quantum maps acting on a simple quantum system. It induces a semi-group 
in the subset of states with maximally mixed partial traces. Subadditivity of the von Neumann 
entropy with respect to this product is proved. It is equivalent to subadditivity of the entropy 
of bistochastic maps with respect to their composition, where the entropy of a map is the 
entropy of the corresponding state under the Jamiolkowski isomorphism. Strong dynamical 
subadditivity of a concatenation of three bistochastic maps is established. Analogous bounds 
for the entropy of a composition are derived for general stochastic maps. In the classical case 
they lead to new bounds for the entropy of a product of two stochastic matrices. 
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1 Introduction 



General quantum dynamics are described by quantum stochastic maps, also called quantum 
channels or quantum operations. It is therefore crucial to investigate their properties in order to 
understand admissible dynamics in a quantum state space. The Jamiolkowski isomorphism [10, 
20] associates with any quantum stochastic operation a quantum state acting on an extended 
space. Up to normalization, this state is nothing else than the dynamical or Choi matrix of the 
map. Thus the features of a quantum map are encoded in a state. 

The spectral decomposition of the Jamiolkowski state yields a canonical Kraus decompo- 
sition of the map, i.e. realizes the map in terms of measurement operators. For a unitary 
evolution, the corresponding state is pure and the Kraus form consists of a single unitary op- 
erator. For the completely depolarizing channel, the state is maximally mixed and the Kraus 



1 



decomposition consists of many terms. The degree of mixing of the measurement operators 
required to construct a quantum map by its canonical Kraus form can therefore be estimated 
by a quantity like the entropy of the map which is actually the entropy of the associated 
state. Thus this entropy vanishes for a unitary evolution and reaches its maximal value for 
the completely depolarizing channel. 

The aim of this paper is to analyze properties of the composition of stochastic maps [22]. 
Using the subadditivity of entropy for composite systems we prove an analogous dynamical 
subadditivity for bistochastic maps, i.e. maps which preserve the identity, see (70). A similar 
inequality proved for the concatenation of three bistochastic maps may be called strong dy- 
namical subadditivity (87). Dynamical subadditivity generalizes to general stochastic maps by 
adding an extra term which vanishes for bistochastic maps. Restricting to diagonal states we 
obtain bounds for the entropy of the product of two classical stochastic matrices, see (106). 
This generalizes a recent result of Slomczyhski on entropy of a product of bistochastic matri- 
ces [17]. 

Composition of quantum maps induces an action in the space of quantum states on a 
bipartite system. We analyze properties of this action which allows to construct a semi- 
group in the space of Hermitian matrices. Composition of states induces also a semi-group 
structure in the set of positive definite operators on a bipartite system whose partial traces 
are proportional to the identity. 

The main tools are coupling techniques, which associate to quantum maps states on com- 
posite systems, combined with subadditivity and strong subadditivity of quantum entropy. 
The explicit constructions that we use are related to similar constructions of Lindblad and 
techniques used in studying quantum dynamical entropy in the sense of [1]. There are also 
connections with quantum coherent information and inequalities as the quantum information 
data processing inequality, see [16]. 

To illustrate the composition of states in action we work out the case of quasi-free Fermionic 
maps. In this setup we derive an explicit form of the composition and discuss dynamical 
subadditivity. 

The paper is organized as follows. In Section II the necessary properties of quantum states 
and quantum maps are reviewed. In particular we consider the case of quantum operations 
with a diagonal dynamical matrix and show that any such stochastic or bistochastic quan- 
tum map reduces to a stochastic or bistochastic matrix acting on classical probability vectors. 
The notion of composition of states is introduced in Section III where some of its properties 
are analyzed. In Section IV we analyze the entropy of maps and formulate dynamical and 
strong dynamical subadditivity for compositions of bistochastic maps. Furtheremore we dis- 
cuss analogous results in a more general case of stochastic maps. More explicit examples of 
low-dimensional quasi-free structures are presented in the Appendix. 
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2 Quantum states and quantum maps 



2.1 Quantum dynamical matrices 

Let p denote a iV-dimensional density matrix i.e. a Hermitian, positive operator, satisfying 
the trace normalization condition Trp = 1. The expectation value of an observable X, i.e. a 
iV-dimensional matrix, is given by the usual relation 

(X)=TrpX. (1) 

Such expectation functionals are called states. We shall in the sequel identify ( ) with its 
corresponding p. Let T>n denote the set of generally mixed quantum states acting on the iV 
dimensional Hilbert space TLn- It is a convex, compact set of real dimensionality N 2 — 1. In 
the case of a qubit, i.e. N = 2, the space of mixed states is the Bloch ball, T>2 = B% C K 3 . 
More generally, there is a one to one correspondence between arbitrary, i.e. not necessarily 
positive linear functionals F on the observables and iV-dimensional matrices p 

F{X)=TrpX. (2) 

General quantum maps, sometimes called super-operators, are linear transformations ei- 
ther of the observables (Heisenberg picture) or of the functionals on the observables (Schrodin- 
ger picture). In this paper we shall mostly use the Schrodinger picture and denote such 
quantum maps by <&. The adjoint & of a super-operator <E> is given by 

*t(p) = ($(p*))*. (3) 

Fixing a basis in Hn we identify a ^-dimensional matrix p with a vector of dimension N 2 
just by writing the entries of p in lexicographical order. So the entry p mil is placed on the 
((m — 1)N + ju)-th row. A general linear quantum map, 

P ^P':=HP) (4) 
may be described by a matrix of size iV 2 still denoted by <£, 

Prnfi = Pnu , (5) 

where Einstein's summation convention is taken. 

Another convenient way to describe quantum maps is to use the Choi-Jamiolkowski encoding 
or dynamical matrix D$ [18]. It amounts to a reordering of matrix elements of $ 

D$ = <I> R SO that (D$)mn = (<E> R )mn = Itnf! . (6) 

7 fxv fxv nv 
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Let us consider the projector on the maximally entangled state |V> + ) := ^7= Y2 m =i \ m ) ® \ m ) 

P + :=\^){i> + \ = ^Y,\ m )^\®\ m )^\ • (7) 



ma 



We can in a similar way as above identify P + with a vector in a space of dimension N 2 x N 2 . 
Its entries are 

^ ] ^fem. $kn $kv = 3mn <^uz/ ■ (8) 
kn 



We now compute 



((*®id)(P + )) = (*®id) (P + ) mV 

$ my 5 nn , 5 VV , (P+) m , n , v 



A m y 

= — $m/i = (D$)mn (9) 
N nv y ' uu 

So we see that, up to a factor AT, is the action of id(£><3? on the one-dimensional projection 
P + on the maximally entangled state. The map $ i-» D$ is linear and it intertwines adjoints 

£>$t = (£>*)*. (10) 

Equation (6) may be considered as the definition of the reshuffling transformation, written 
<3? — > <I? R , which is defined for any matrix $ acting on the Hilbert space 1~i^2 = TIn ®TIn [20]. 
It should be stressed that the reshuffling operation depends on the distinguished basis in TCn 
that we have used. 

Choi's theorem [4], proves that a map <!> is completely positive (CP), which means that the 
extended map $ ® id is positive for any size of the extension, if and only if the corresponding 
dynamical matrix D<j>, also called Choi matrix, is positive, D$ > 0. The eigenvalue decompo- 
sition of the dynamical matrix of a CP map <3? leads to the canonical Kraus form [11] of the 
map 

N 2 

®{p) = Y,A aP Al. (11) 

a=l 
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where the Kraus operators are orthogonal 

(A a \Af,) = TtA$ x A l3 = d a 8 a { ) , (12) 

so that the non-negative weights d a become the eigenvalues of the dynamical matrix 
Hence, in this (almost) canonical form the number of Kraus operators does not exceed iV 2 . 

A quantum map $ is trace preserving (TP) if Tr$(p) = Trp for any p. The corresponding 
dynamical matrix acts on the composite Hilbert space H. N 2 = Ha ®Hb and, in terms of 
the dynamical matrix, trace preserving means that 

Tr A D* = l. (13) 

This implies in particular that TrD$ = N. Completely positive trace preserving maps (CPTP 
maps) are often called quantum operations or quantum stochastic maps. Since the dynamical 
map of a quantum operation $ is positive and normalized as in (13) the rescaled matrix Z?$ 
is a state on the extended Hilbert space Hn <8> 7~Ln, see [10, 3]. We shall say that 

s--=jrD* (14) 

is the Jamiolkowski state associated to <&. A quantum map is called unital if it leaves the 
maximally mixed state invariant. This is the case iff 

Tv B D^ = t, (15) 

a condition dual to (13). A CP quantum map which is trace preserving and unital is called 
bistochastic. A composition of bistochastic maps is still bistochastic. 



2.2 Classical case - diagonal dynamical matrices. 

Let us diagonalize a density matrix. The elements on the diagonal may be interpreted as a 
classical probability vector P with components pi = pa, i.e. pi > and YliLiPi = 1- In a 
similar way a quantum map reduces to a classical one if the dynamical matrix D = <& R is 
diagonal, D ab = T ab 5 ac 5 bd . 

cd 

Reshaping the diagonal of the dynamical matrix which has dimension A^ 2 one obtains a 
matrix T of dimension N 

T = T{$) with T ij = ^ ii , (16) 

33 

with no summation performed. Positivity of D implies that all elements of T are non negative. 
Furthermore, the partial trace condition (13) implies that the matrix T is stochastic since 
YliLi Tij = 1 f° r all j = 1, 2, . . . , N. In fact, on diagonal matrices the action of a diagonal 
dynamical matrix D reduces to a Markov transition of a probability vector, P' = TP. 
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If, additionally, the complementary partial trace condition (15) holds, then the matrix T 
is bistochastic Ylf=i ^ij = 1 for all i = 1, 2, . . . , JV and the uniform vector P* := (-^, . . . , -^) 
is invariant under multiplication by a bistochastic matrix T. Hence quantum stochastic and 
bistochastic maps acting in the space T>n of quantum states can be considered as a direct gen- 
eralizations of stochastic and bistochastic matrices, which act on classical probability vectors. 

Note that the quantum identity map, <I> = id, is not classical, since the dynamical matrix 
D = <£ R is not diagonal, and the off-diagonal elements of p are preserved. On the other hand, 
the coarse graining map, 

N 

*CG(p):=X>><t|p|i><i| (17) 
i=i 

which strips away all off diagonal elements of a density operator is described by a diagonal 
dynamical matrix. The corresponding stochastic matrix is the identity, T(<I>cg) = 1, since all 
diagonal elements of p remain untouched under the action of (17). Let us also distinguish the 
flat stochastic matrix T* whose elements are all equal, (P*)ij = jj- It is a bistochastic matrix 
which maps any probability vector P into the uniform one, T*P = P*. 

2.3 Quasi-free Fermionic states and maps 

Quasi-free states and maps on a Fermionic algebra will be used as an example. More details 
and references to the original papers can be found in [1] and [7]. In appendix we shall describe 
in more detail these objects for systems with few modes. 

The algebra A(Hn) of observables of a iV-mode Fermionic system is generated by creation 
and annihilation operators a*(ip) and a(ip) where ip belongs to the one-particle space Hn- 
These operators satisfy the canonical anti-commutation relations (CAR) 

a*(tf+atp) = a*((p)+aa*(ip) , {a* (tp) , a* (ip)} = and {a*(<p), a(ip)} = (ip,(p) ■ (18) 

It is not hard to see that A{T~Ln) is isomorphic to the algebra of matrices of dimension 2 N . 
The space Hat is called the one-particle space. 

There exists a widely used class of states and maps called quasi-free. These objects corre- 
spond quite literally to Gaussian objects for Fermionic systems and they are fully characterized 
by operators on the one-particle space. Essentially every effective description of Fermionic sys- 
tems, such as the Hartree-Fock approximation, is based on quasi-free structures. 

Proposition 1 (Quasi- free states). Let Q be a N -dimensional matrix. Every element of 
A(H.n) can be written as a linear combination of ordered monomials, i.e. monomials in creation 
and annihilation operators where the a* 's appear to the left of the a 's. Define a linear functional 
( )q on A(TLn) by extending linearly its definition on ordered monomials 

(a%<p 1 )---a*(<p n )a(il> n )---a(il> 1 )) Q := Det ( [(^ , Q <fi)] ) (19) 
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and 



(mon)g = for every other ordered monomial mon. (20) 

Then ()q is a state on A(Hn) iff < Q < 1, i.e. Q is Hermitian and every eigenvalue of Q 
belongs to (0, 1). States of the type ( )q are called quasi-free and Q is called the symbol of the 
state. 

The set of quasi-free states on A(TCn) is not convex, so it is natural to consider the convex 
hull of the quasi-free states which is a proper subset of the full state space. It is known that 
any quasi-free state can be decomposed into a convex combination of quasi-free states whose 
symbols are projectors. Moreover, a quasi-free state is pure iff its symbol is a projection 
operator. Therefore the extreme points of the convex hull of quasi-free states consist of the 
states with a projector as symbol. It is remarkable that the extreme points of the convex set 
of symbols (0,1) precisely consists in the projection operators on Hn- The maximally mixed 
state, i.e. the normalized trace on A(Hn) is quasi-free and given by Q = \ 1. 

Let A(T~Cn © 7~Cn) be a CAR algebra for a composite system, it is isomorphic to the tensor 
product of two copies of A(TLn). It is however, because of the Fermionic nature of the system, 
more natural to consider the graded tensor product where each of the two parties is identified 
by the maps 

a*(v?i) ^ a*(</?i © 0) and a* ((p 2 ) ^ a* (0 ® (p 2 ) ■ (21) 
The quasi-free state with symbol 

1 (I 1) <*> 

is pure on A{TLn © ~Hn) and its restriction to each of the parties A(Hn) is maximally mixed, 
therefore it is a maximally entangled state on the composite system. 

Quasi-free CPTP maps are conveniently described in Heisenberg picture by their explicit 
action on monomials [6, 9]. For the purpose of this paper it will suffice to know their action 
on quasi-free states. 

Definition 1 (Quasi- free CPTP maps). A CPTP quasi- free map is a CPTP map on the state 
space of A(TCn) which maps quasi-free states in quasi-free states. 

It can be shown [5] that quasi-free TPCP maps are determined by a couple of iV-dimensional 
matrices R and Z which satisfy the constraint 

< Z < 1 - R*R . (23) 

Moreover, the action on quasi-free states is given by 

()q^( )r*qr+z ■ (24) 
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It is easily seen that the constraints (23) are necessary and sufficient to guarantee that for 
every symbol Q the expression R*Q R + Z is again a symbol. In order to mark the dependence 
of the map on R and Z we shall use the notation &r^z- The map ( )q i— > ( )r*qr+z is affine 
on (0, 1) but not linear but we may use the standard representation 



Q 



and ^-(T'fd)' (25) 



Here Ad(i?) is the adjoint map X R*X R. Composing $>r 2 ,z 2 ° ^R 1 ,z 1 returns the map 
®R 3 ,Z 3 with 

R 3 = R 2 R X and Z 3 = R\Z X R 2 + Z 2 . (26) 
In terms of the representation (25) this amounts to the usual matrix multiplication 

/Ad(i? 3 ) Z 3 \ _ (Ad(R 2 ) Z 2 \ /Ad(fli) ZA 

V o id) ~ V o id J V o id) ■ 120 

The representation (25) is the quasi-free analogue of (5). 

Using (22) we see that the Jamiolkowski state associated with $r 5 z is quasi-free on A{Hn® 
TLn) with symbol 



l 



1 R 
R* R*R + 2Z 



(28) 



The set of quasi-free CPTP maps is again not convex. A quasi-free map &r : z is extreme 
within the set of CPTP maps iff 



Z = Vt- R*R Py/1- R*R (29) 

with P a projection operator. The same construction that applies for quasi-free states allows 
to decompose any quasi-free CPTP map into a mixture of extreme quasi-free CPTP maps. 
Therefore the extreme points of the convex hull of quasi-free CPTP maps consists precisely of 
the maps of the form (29). 

Finally, a quasi-free CPTP map <&R t z is bistochastic iff 

Z = \ (1 - R*R) . (30) 

This follows from the invariance of the tracial state which has symbol \ 1. As such a map is 
fully determined by a single A^-dimensional matrix R with ||J?|| < 1, we simply write <&r. 

2.4 Entropies of maps and states 

To any normalized probability vector P of size we may associate its Shannon entropy 

N N 

H(P) :=-J>kiPi = I>(P0 , (31) 

i=i i=i 
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where we have introduced the function 

rj(x) := — xlnx for x > and n(0) : = . (32) 

This entropy is a measure for the mixedness of a probability vector. In a similar way the 
degree of mixing of a quantum state p is characterized by its von Neumann entropy 

S(p) :=-Trpmp = Trr/(p) , (33) 

equal to the Shannon entropy of its spectrum. The entropy varies from zero for a pure state 
to In TV for the maximally mixed state, p* = 1. 

The density matrix <j = D^/N associated to a quantum stochastic map <3? depends on the 
basis that is used to compute the entries of see (5). A change of basis in Hn corresponds 
to a unitary transformation of -D$, therefore the eigenvalues of D§ don't change. It is then 
natural to consider the von Neumann entropy of the bipartite state associated with <J> by the 
Jamiolkowski isomorphism. 

Definition 2. Let be a trace-preserving completely positive map with associated Jamiolkow- 
ski state q = D<&. The entropy of the map $ is defined to be 

S(*):=S(?) = S(^D*) . (34) 

Since the state ? = D$/N acts on the extended Hilbert space TLn^TLn the entropy of the 
map varies from zero for a unitary dynamics to 21niV for a completely depolarizing channel 

which sends any state to the maximally mixed state, $*(p) = p*, see [20]. 

Let us now move to a classical discrete dynamics in the probability simplex. The following 
definition of entropy of a stochastic matrix introduced in [17, 21] 

N 

H!(T):= J>Jh%) , (35) 
j'=i 

is decorated by a label "I", as it is based on an invariant state of a matrix, P l = TP 1 . Here 
tj denotes the j-th column of a transition matrix T, so (35) represents the average Shannon 
entropy of columns of T weighted by its invariant state P 1 . 

To demonstrate a direct relation to the quantum dynamics we shall use a simplified version 
of the entropy of a transition matrix 

N N 

HCO^-tfEE 2 * 1 * 1 *- (36) 

i=l 3=1 

Observe that for any bistochastic matrix the uniform vector is invariant, P 1 = P*, so = jj 
and both definitions of entropy do coincide. Both quantities, Hi(T) and H(T), vary from zero 
to In N. 
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Using eq. (16) one concludes that for any stochastic map <I> represented by a diagonal 
dynamical matrix D<j> its entropy is up to a constant equal to the entropy of the associated 
stochastic matrix, 

S(<E>) = H(T($>)) + In AT . (37) 

The constant In TV is due to the normalization factor in front of the dynamical matrix, see 
definition (34). Although in general the entropy of a quantum map belongs to [0, 2 In AT], the 
entropy of the maps represented by diagonal D and corresponding to stochastic matrices vary 
from lniV to 21niV. Among this class the minimal entropy characterizes the coarse graining 
map (17), for which T(<J>cg) = 1> and S( < I ) cg) = hiiV. The maximum is achieved for the 
completely depolarizing channel, since T(®*) = T* so that S(3>*) = H(T*) + ln iV = 21n N. 

The von Neumann entropy of a Fermionic quasi-free state with symbol Q can easily be 
expressed in terms of Q 

S qf (Q) = Tr(r / (Q)+r / (l-Q)) . (38) 

Quite explicit expressions can be given for the entropies of extreme CPTP and bistochastic 
quasi- free maps using (28) and (29). 

For an extreme CPTP quasi- free map, we obtain after some algebraic manipulations 

S(<t> R , z ) = S^(l(t + \R\ 2 -2\R\P\R\)) 

= TV^(i(l + |i?| 2 - 2\R\ P \R\)) + r/(i(l - |i?| 2 + 2\R\ P |i?|))) . (39) 

In these expressions P is a projector and Z = \Jt_ — |i?| 2 P y/l — \R\ 2 - 
For a bistochastic map determined by R with ||i?|| < 1, we obtain 

S(* fl ) =2^(^(1(1 + \R\)) + r?Q(l - \R\))) . (40) 
2.5 Entropy exchange and Lindblad's theorem 

Consider a CP map $ represented in its canonical Kraus form (11). For any state p 6 D/v 
define a positive operator a = &(&,p) acting on the extended Hilbert space TCn2, 

d afi :=Tr pA^A a , a, (3 = 1, . . . , N 2 . (41) 

If the map $ is stochastic, then the operator a is also normalized in the sense that 

N 2 

Tr a = TrpAlA a = Tr p = 1 (42) 

a=l 

and so it represents a density operator in its own right, a G T> N 2. In particular, if p = p* = 1 
then, using the canonical Kraus decomposition (11), one shows that 

<r(cD, /0 ,) = i-^= ? . (43) 
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The von Neumann entropy of a depends on p, and equals S(<5), as denned above, if p is the 
maximally mixed state. 

Auxiliary states a in an extended Hilbert space were used by Lindblad to derive bounds 
for the entropy of the image p' = <&(p) of an initial state under the action of a CPTP map. 
Lindblad's bounds [13] 

|S(<T)-S(p)|<S(p')<S(a) + S(p), (44) 
are obtained by defining yet another density matrix in the composite Hilbert space TLn <S> 

M M 

(45) 

a=l 0=1 

where M = N 2 and {\a}} is an orthonormal basis in Hm- Computing partial traces one finds 
that 

Tr/v w = a and Tr m w = P ■ (46) 

It is possible to verify that S(u>) = S(p), and so one arrives at (44) using subadditivity of the 
entropy and the triangle inequality [2]. These results can be obtained using the first part of 
the proof of Theorem 1. 

If the initial state is pure, that is if S(p) = 0, we find that the final state p' has entropy S(<r). 
For this reason S(<r) was called the entropy exchange of the operation $ by Shumacher [15]. In 
that work an alternative representation of the entropy exchange was given 

S(a(<M) =S((id®<&)| ¥ >>< ¥ >|) , (47) 

where \ip) is an arbitrary purification of the mixed state, Trp \<p)(ip\ = p. To prove this useful 
relation it is enough to find a pure state in an extended Hilbert space, such that one of its 
partial traces gives a and the other one the argument of the entropy function at the right 
hand side of (47). 

A kind of classical analogue of the quantum entropy bound (44) of Lindblad was proved 
later by Slomczyhski [17]. He introduced the notion of entropy of a stochastic matrix T with 
respect to some fixed stationary probability distribution P = {pi}f =1 , 

N 

H P (T) := Y^Pi H (^) £ = (T u ,T 2i , . . . , T Ni ) . (48) 
i=i 

This quantity — an average entropy of columns of matrix T weighted by probability vector 
P — allows to obtain the bounds for the entropy of a classically transformed state, P' = TP, 

Hp(T) < H(P') < H P (T) + H(P) . (49) 
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These bounds look somewhat similar to the quantum result (44) of Lindblad, but a careful 
comparison is required. Applying the definition (41) of the auxiliary state a to the classical 
case of a diagonal state, pu = Pi&ij and a diagonal dynamical matrix D we find 

S(a class ) = H P (T) + H(P) , (50) 

where T = T(D R ) is the classical stochastic matrix given by (16). Substituting this result into 
the Lindblad bound (44), and renaming p and p' into P and P' we realize that the argument 
of the absolute value in the lower bound reduces to Hp(T) and is not negative, so we arrive at 

Hp(T) < H(P') < Hp(T) + 2H(P) . (51) 

The lower bound coincides exactly with the result (49) of Slomczyhski. The upper bound is 
weaker (note the presence of the term 2H(P) instead of H(P)), but it holds in general for all 
quantum maps, while (49) is true for classical dynamics only. 

3 Composition of maps and composition of states 

We first recall some properties of the reshuffling transformation of a matrix as defined in (6). 
Reshuffling does not preserve the spectrum nor the Hermiticity of a matrix. It is an involution, 
since performing this transformation twice returns the initial matrix, (X R ) R = X. 

Using the Jamiolkowski isomorphism the composition of maps acting on a single quantum 
system can be used to define a composition between quantum states of a bi-partite system, in 
fact, this composition extends to arbitrary matrices. 

Definition 3. The reshuffling operation (6) defines a composition between arbitrary matrices 
o\ and (72 on a composite system TLn ®~Hn 

ai Qa 2 := K(J 2 R ) R . (52) 

For stochastic matrices obtained by reshaping two diagonal density matrices T\ = T(a R ) 
and T 2 = T(a R ) according to (16), the composition of the diagonal states returns the usual 
multiplication of stochastic matrices, 

T((<7i a 2 f) = T(a R ) T(o~ 2 ) = T X T 2 . (53) 

Using the definition of reshuffling we see that generally (D R ) 2 differs from (D 2 ) R . Therefore 
the composition performed on two copies of a state a differs from its square 

a® 2 ■= a a + a 1 . (54) 



12 



3.1 Properties of the composition 

Lemma 1. Let X and Y denote two matrices of size N 2 . 

IfX>0andY>0 then X Y = {X R Y R f > . (55) 

Proof. Denoting by ^ and <I> the completely positive maps with corresponding dynamical 
matrices X and Y we see that X Y is the dynamical matrix of the composed map ^ o 
Since the composition of two completely positive maps yields again a completely positive 
map [3], we infer (55). □ 

Proposition 2. The set of all operators acting on a composite Hilbert space T~Ln 1~Ln > 
equipped with the composition law (0), is a non-Abelian associative semi-group. Moreover, if 
o\ and o"2 are Hermitian operators on TL^ ®~Hn then also o\ 02 is Hermitian. Therefore 
the set of all Hermitian operators on TLn Ti-N is a non-Abelian associative sub semi- group. 

Proof. Let P + denote the projector on the maximally entangled state. Since (P+) R = 1, this 
operator plays the role of the neutral element of the composition. The composition is non- 
Abelian, o~\ 0d2 7^ o~2 ffi, because the composition of quantum maps is not commutative. 
It is on the other hand associative (a± 02) 03 = o\ (02 03), since the composition of 
maps is. 

It remains to prove is that o~\ 002 is Hermitian if o~\ and 02 are. Let $ be a super-operator 
associated with D<j> through the Jamiolkowski isomorphism, then = D^. Moreover, 

for two super-operators <I> and ^ we have 

(*o*)t(X) = (<f>oV(X*))* = = &((y{X*))*) = (&otf)(X) (56) 

Using 

£>$o# = D$ (57) 

finishes the proof. □ 

Restricting our attention to the set V N 2 of quantum states on a composite system, we see 
that this algebraic structure breaks down since the trace condition, Trcr = 1, is not preserved 
under composition. However, one may overcome this difficulty by selecting a certain subset 
of quantum states. Thus consider the subset of density matrices of a composite system 
Ha §Wb of size N 2 such that their partial trace over the first system is the maximally mixed 
state 

V\ 2 := {a e V N 2 : Tr A a = ± 1} . (58) 

With respect to the Jamiolkowski isomorphism these states correspond to the trace preserving 
maps. Since a composition of any two trace preserving maps preserves the trace we infer that 
the composition (52) acts internally in the set V l N2 . 
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Proposition 3. The set of all operators acting on the Hilbert space Hn of a bipartite system 
such that their left marginal is proportional to the identity, equipped with the composition 
is a non-Abelian associative semi-group. 

It is convenient to distinguish another composition sub-algebra by defining the set of states 
with both marginals proportional to the identity, 

V% 2 := {a e V l N2 : Tx B a = ± 1} . (59) 

Due to condition (15) this semi-group is generated by compositions of bistochastic maps. 

3.2 Idempotent states 

Consider the state a of a bipartite system obtained by extending an arbitrary state p by the 
maximally mixed state 

<r:=^p®l. (60) 

This state is proportional to the dynamical matrix D of the operation $ p , which acts as a 
complete single-step contraction, sending any initial state to into p, 

a = jj D<$, p , where $ p {uj) := p (61) 

for any w e To show this let us start with the dynamical matrix of this map, Dmn = 
Pm/j, o~nv ■ Writing out the matrix entries of 

J := & p (u) = D R uj = (p®±) R uj (62) 

in the standard basis we obtain the desired result 

u'mti = DmnLO nv = p mM (Tr Uj) = p mfJi . (63) 

A state a of an algebra is called idempotent if aQo = a. Hence any state of the form (60) is 
idempotent, since the single-step contracting map Q p applied a second time does not influence 
the system anymore, <3? p o $ p = $ p . 

3.3 Composing quasi- free quantum maps 

The Jamiolkowski state associated to a quasi-free quantum map <&r,z has the symbol 

°^(^ M + 2z)^ 1 - (64) 
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When R and Z vary through the matrices of dimension N subject to the constraint < Z < 
1 — R*R the symbol (64) varies trough the matrices of dimension 2N in the interval (0,1) 
which are of the form 

» U S) • (65) 

Remark that this set is convex. Composing quasi-free quantum maps induces now the following 
composition law on these matrices 

1 (t R 1 \ , / 1 SA 1 ( l rA 

2 v^i 2 V 5 * -v ' 2 V r r t 2 ) 

with 



Ti = 12i5i and T 2 = S^ifc - + S 2 . (66) 

This composition law is the analogue of (52). It is associative and non-commutative but only 
affine in its first argument. 



4 Entropy of a composition 

In this section we analyze the behaviour of the entropy of a quantum operation under com- 
position. The bounds (49) on the increase of entropy of probability vectors under discrete 
dynamics allowed Slomczyhski to prove the subadditivity relation [17, 21] 

Hi(Ti) < Hi(T 2 7i) < Hi(Ti) + Hi(T 2 ) (67) 

provided both stochastic matrices T\ and T 2 have the same invariant state, P\ = P\. Restrict- 
ing our attention to the case of bistochastic matrices for which P\ = P\ = P* = • • • , ^} 
we use (36) instead of (35) and drop the label "I" in the subadditivity relation to get 

H(Ti) < H(T 2 Ti) < H(Ti) + H(T 2 ) . (68) 

Considering a product of three bistochastic matrices Slomczyhski proved [17] a strong subad- 
ditivity relation for classical dynamics, 

H(r 3 r 2 ri) + h(t 2 ) < h(t 3 t 2 ) + h(t 2 Ti) . (69) 

4.1 Dynamical subadditivity and strong subadditivity for bistochastic maps 

Motivated by the classical results above we formulate and prove their quantum counterparts. 
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Theorem 1 (Dynamical subadditivity for bistochastic quantum operations). Let <I>i be a 

bistochastic quantum operation and $2 a general stochastic quantum map then their entropies 
satisfy the subadditivity inequality 

S($ 2 °$i) < S($i) +S($ 2 ) • (70) 

If both $1 and $2 are bistochastic then 

max({S($i), S($ 2 )}) < min({S($i o $ 2 ), S($ 2 ° $1)}) • (71) 

-An equivalent statement of (71) is the triangle inequality for composition 

max({S(<7i),S(<7 2 )}) <mm({S(a 2 Qa 1 ),S(a 1 Qa 2 )}) < S(<7i) + S(<7 2 ) , (72) 

where o-\-,o~2 G and the set has been defined in (59). 

Proof, i) We first show the upper bound (70). 

Let be a stochastic quantum operation, generally not bistochastic, with Kraus form 

K 

e(p) = Y,C aP Cl. (73) 

a=l 

Introduce a map from TCn to TCn <8> TCn ® as follows: fix an orthonormal basis {\a)} in 
TCk and let 

K 

F V :=Y,{C aV )®\a) , (74) 

The adjoint map acts as 

FV ® |q) = (75) 

and one checks that 

K K 

F^Fp = Ft J2( C ^) ® 1$ = £ = V > (76) 

0=1 0=1 

since is trace preserving. Therefore F is an isometry. It follows in particular that, for an 
arbitrary iV-dimensional matrix p, F pF^ and p have up to multiplicities of zero the same 
eigenvalues. 

Using F we can express the Lindblad operator ui, see (45), as 

K 

cu = FpFt = C aP Cl \a)(P\ . (77) 

Q,/3=l 
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The operator u is a density matrix on the composite system Hn <8> T~Lk- If the initial state 
is maximally mixed, p = p*, taking the partial trace over the first subsystem we obtain the 
following density matrix of the ancilla E = Hk , 

p = TrNUJ = a (0 :P ^ = ? = 1 D@ . ( 78 ) 

We perform the construction of above for the composition $2 We first write both 
quantum operations $1 and $2 in their Kraus forms 

M M 

$ 1 (p) = ^A Q pAt and $ 2 (p)=J]5 Q p4 (79) 

Q=l Q=l 

where M = N 2 . Putting = $2 $1 we consider its Kraus decomposition 

C Q2Ql := B a2 A ai (80) 

implying that K = N 2 x N 2 = N 4 . Consider a state U21 acting on an extended tripartite 
space Hn <8> % <8> Hm = Hn ® E 2 ® E\ 

N 

0J21 = ^2 C a2ai p C^ 2/3i ® |a 2 ®ai)(/?2®/3i| 

oi,02,/3i,/32=l 
AT 

® |a 2 ®ai)(/3 2 ®/3i| (81) 

"i 5 Q!2,/3i,/32=l 

and let P21 be the restriction of uj 2 \ to the ancilla E 2 E\ = Hm <8> %■ We also compute the 
restrictions pi and p2 to the first and second ancilla. Assuming that p = p* = -^t we obtain: 

P12 = ^$20$! (82) 



Pi = Tr nN ® E2 U21 = Tr nN — ^ A^tA^ ®\a 1 ){p 1 \ = — (83) 

Ql,/9l = l 

1 ^ 1 

P2 = Tr WiV0Bl = Tr WiV — B^lB^ ® |a 2 )(/? 2 | = — D$ 2 . (84) 

C«2,/?2=l 

The last inequality holds because 3>i is bistochastic. The upper bound of (70) now follows by 
applying subadditivity of the entropy to the state P12, see [14, 3]. 

ii) The lower bound is a special case of the bound in Theorem 3. As $1 and $2 are 
bistochastic 

$i(p*) = $ 2 (p*) = P* (85) 

and the contribution of the terms within the square brackets vanishes. □ 
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Corollary 1. Let & be a bistochastic map, then 

S($ on ) < raS($) . (86) 

It should be remarked that (86) is only meaningful when S(<1>) <C IniV and n is not too 
large as the entropy of a map is anyway bounded by 2 In TV. 

In a similar way one may analyze properties of a concatenation of three consecutive oper- 
ations. Motivated by the inequality (69) for the entropy of the products of three bistochastic 
matrices, we formulate its quantum mechanical counterpart. 

Theorem 2 (Strong dynamical subadditivity for bistochastic quantum operations). Let $i ; 
<E> 2 and $3 be bistochastic quantum operations then their entropies satisfy the inequality 

S($ 3 o $ 2 o $i) + S($ 2 ) < S($ 3 o $ 2 ) + S($ 2 o $i) . (87) 

This is equivalent to 

S{a 3 0-2 <7i) + S(<7 2 ) < S(cr 3 <r 2 ) + S(a 2 a x ) , (88) 
w/iere ai,a 2 ,cr 3 £ V^ 2 . 

Proof. Consider three bistochastic quantum operations <J>i, <& 2 and <I>3 acting in sequence on 
the maximally mixed state p* £ Mn- We repeat the construction of the Linblad operator as 
in (77) but now for three ancillas Ei, E 2 and £3 and obtain a density matrix acting on 
Hn -E3 -E 2 E\. The restrictions of W321 to some of the ancillas will be denoted as before 
by ^'s. We now write the strong subadditivity of quantum entropy for the system E^E 2 E\ [12] 



S(p 3 2l) + S(p 2 ) < S(/9 2 l) + S(/> 3 2) • (89) 

This yields, using the bistochasticity of $3 and <3?2 

S(p 2 )=S($ 2 ) (90) 

S(p2i) =S($ 2 o$i) (91) 

S(p 32 ) =S($ 3 °$ 2 ) (92) 

S(p 3 2l) = S($ 3 °$2°$l) , (93) 

which ends the proof. □ 



4.2 Generalization of dynamical subadditivity for stochastic maps 

Results obtained in previous section for bistochastic maps can be generalize for the case of 
arbitrary stochastic maps. 
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Theorem 3 (Dynamical subadditivity for quantum operations). Let $1 and $2 be quantum 
operations then their entropies satisfy the inequalities 

S(*i) + Ai < S($ 2 o $1) < S(*i) + S($ 2 ) + A 2 (94) 

where 

Ai = S($ 2 o - S(*i(p*)) , (95) 

A 2 = s(a(* 2 ,^i(^)))-S(* 2 ). (96) 

Proof. We repeat the proof of Theorem 1 for the composition $2 $1 up to the construction 
of the state uj on TCn ® £2 <8> £1 , see (77). If the first operation is not bistochastic then the 
restriction (84) has the form 

N 

P2 = Tt Hn ® Ei L021 = Tt Hn 5a 2 *i(p*)^® |a 2 )<#2| = <t($ 2 ,$i(p*)) ■ (97) 

012,^2=1 

The restrictions P12 and pi are the same as before. The entropy of a forms the first term of 
A2 so subadditivity of the entropy of P12 implies the upper bound of (94). 

To prove the lower bound consider (81). As an initial state let us take p* = ^1 and purify it 
by adding an additional party R of dimension TV. This yields a pure state on R®Hn®Ei®Ei 
defined by the normalized vector 

j N N 2 

nYI J2 I*)® ( B «2^il*)) ®l<*2>® • (98) 

i=l ai, 012=1 

We now use the strong subadditivity of the entropy denoting by pr the restriction of this pure 
state to the party R with similar notations for restrictions to other parties 

S{pre 2 e 1 ) + S(p El ) < S{p REl ) + S{pe 2 e 1 ) • (99) 

Computing all the terms appearing in (99), 

S(p El ) =S(*i) (100) 

S(pE 2El ) = S($ 2 °$l) (101) 

S( PREl ) = S^M) (102) 

S( P RE 2El ) =S($2 0$lW) (103) 

we arrive at the lower bound (94). □ 
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The above inequalities formulated in the language of quantum maps obviously hold for the 
composition of states which belong to the set (58) and its subset (59). 

The above results may be linked to properties of coherent information defined by Schu- 
macher and Nielsen [16] as a function of exchange entropy S(cr(<]?,p)): 

I($,p) = S(*(p))-S(a(*,p)). (104) 
Denoting by I\ the coherent information for the first operation, I\ = S(3>i(p)) — S(<7(<]?i, p)) , 
and by hi the analogous quantity for the concatenation hi = S ^3>2(^i(p))^ — S(&($>2°&1, p)) , 
these authors proved [16] the quantum data processing inequality: 

hi < h ■ (105) 
This result implies the lower bound of (94). 

4.3 Entropy of a product of stochastic matrices 

Restricting our attention to the case of diagonal dynamical matrices we may analyze classical 
analogues of the above results. In this case an application of a quantum map corresponds to 
action of a stochastic matrix T on a classical probability vector P. 

Theorem 3 implies that in the classical case for an arbitrary stochastic matrices T\ and T 2 
the following bounds hold, 

H(Ti) + [h(T 2 Ti(P,)) - H(Ti(P,))] < H(T 2 Tx) < H(Ti) + H TlP ,(T 2 ) + H(TiP,) , (106) 

where P, = £(1,...,1). 

However it is possible to obtain a stronger upper bound by using properties of strong sub- 
additivity of entropy. Inequality (89) is equivalent to an inequality for the exchange entropy, 

S(a($ 3 O*2O$l,/0*))+S(ff($2,*l(p*))) < S(<7($ 2 0$l,p*))+S(<7($ 3 0$2,$l(p*))) • (107) 

By restriction to the classical case we obtain an inequality for three stochastic matrices T\ = 
7i($i), T a = T a (<I> 2 ) and T 2 = ^($3) (note the notation chosen, which is convenient to state 
the final result). This inequality is similar to the strong subadditivity, but different weights 
for entropies are used. By substitution T a = 1 we arive with a result analogous to (94). 

Theorem 4. Let T\ and T 2 be arbitrary stochastic matrices. Then the entropy of their product 
is bounded by 

H(Ti) + <Ji < H(T 2 Ti) < H(T 2 ) + H(Ti) + S 2 (108) 

where 

< 5i = H(T 2 T 1 (p t ))-H(T 1 (p t )), (109) 
<5 2 = H TlPt (T 2 )-H(T 2 ) . (110) 
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This classical version of (94) valid for arbitrary stochastic matrices T\ and T 2 can be 
considered as a direct generalization of the result of Slomczyhski [17] obtained for bistochastic 
matrices. If both matrices T\ and T 2 are bistochastic a complementary lower bound for the 
entropy of their product holds H(T2) < H(T2Ti), so in this case (68) can be rewritten in the 
stronger symmetric form [21] 

max({H(Ti), H(T 2 )}) < min({H(r 1 T 2 ), H(T 2 Ti)}) . (Ill) 

4.4 Dynamical subadditivity of quasi-free bistochastic quantum operations 

The general results of Theorems 1 and 3 can of course be applied to the case of quasi-free 
completely positive maps. It is however also possible to derive them directly for such special 
maps bypassing e.g. the strong subadditivity that was used in the general proof. By way of 
illustration we prove Theorem 1 for this class of maps. 

Let R\ and i? 2 be iV-dimensional matrices such that ||i?i|| < 1 and \\R2W < 1- Remark 
first that 

<5>R 2 = ®R 2 Ri ■ ( 112 ) 

As 

= R2R1R1R2 ^ R2R2 = (113) 

and as 

x G [0, 1] 1 ^ r,{\ (1 + x)) + r,{\ (1 - x)) (114) 
is monotonically decreasing, we obtain 

S(cD fl2 ) < S($ R2Rl ) = S{$ R2 o $ Rl ) . (115) 

Using S(<£ft) = S(<E>j?*) allows to exchange the roles of R\ and i? 2 and so 

max({S(<I> Rl ),S(<I> R2 )}) <S(^ 2 o$ Rl ) . (116) 

If we denote by {\i(R), A 2 (i?), . . . , Xn{R)} the singular values of the TV-dimensional matrix 
R arranged in decreasing order, then 

Ai(i? 2 i?i) 

> max^A^) Ajv- j+ i(i?i) : 1 < j < AT}) 

> max(^{X jl (R2) Ajv-^+iORi) + A i2 (i? 2 ) Ajv-^+i^i) : 1 < ji < h < JV}) 

i (117) 



21 



It follows that 

{\ (1 + \j{BaRi)) : 1 < j < N} \J{± (1 - Xj(R 2 Ri)) : I < j < N] 
is less mixed than 

{i(l + A J (i? 2 )A JV _ j+ i(i?i)) : l<i< N}(jq(l-\ j (R 2 )X N - j+ i(R 1 )) : 1 < j < N] . 

(118) 

This implies by concavity of r\ that 

N 

S($ R2Rl ) < 2^(r?(i (1 + Aj(-R 2 ) Ajv-j+iCiJi))) + r/(± (1 - A,(i? 2 ) A JV _, +1 ( J R 1 )))) . (119) 

3=1 

Finally observing that for < a, b < 1 

{|(l + o6),|(l-a6),0,o} (120) 

is less mixed than 

{ \ (1 + a), \ (1 - a)} x [\ (1 + 6), i (1 - b)} (121) 

we obtain 

S($ R2 o<l> Rl )<S($ Rl )+S($ R2 ) . (122) 
5 Concluding remarks 

In this paper we introduced the composition between states on a bipartite system and 
derived some basic properties. This composition reflects the concatenation of quantum oper- 
ations under the Jamiolkowski isomorphism. Next, we introduced a simple notion of entropy 
of a quantum map in order to quantify its randomizing properties. We proved the property 
of subadditivity for a composition of arbitrary bistochastic maps and found its generalization 
for the case of stochastic quantum maps. The connection between maps and states allows to 
formulate these properties purely in terms of states and the composition. A restriction to 
the classical setting leads to a generalization of recently obtained bounds on the entropy of 
a product of two bistochastic matrices [17] to the case of a product of arbitrary stochastic 
matrices. 

Recently [19], the concatenation of quantum maps has been investigated from the point of 
view of divisibility. The authors consider the determinant of a super-operator instead of the 
entropy and show that it is contractive with respect to composition. 

A more detailed understanding of the randomizing properties of a quantum map should 
be provided by constructing a Markov like process. A possible track is to generate a state on 
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a spin half-chain in the spirit of finitely correlated states [8] and to consider the associated 
dynamical entropy. This is a subject of future research. 

Acknowledgements: We enjoyed fruitful discussions with I. Bengtsson, V. Cappellini and 
W. Slomczyhski. We acknowledge financial support by the bilateral project BIL05/11 between 
Flanders and Poland, by the Polish Ministry of Science and Information Technology under 
the grant 1 P03B 042 26 and by the European Research Project SCALA. 



A Appendix 

The aim of this appendix is to make for low dimensional one-particle spaces the connection 
between quasi-free states the usual density matrix description of a state more explicit. The 
CAR algebra An '■= A(Hn) is isomorphic to the algebra of matrices of dimension 2^. Be- 
cause of the finite dimensionality all representations are equivalent. A particularly useful 
representation is on Fermionic Fock space 

t(h n ) :=cen N en N ^ s ®-®c . (123) 

Here H N ' as denotes the subspace of antisymmetric vectors under the action of the permutation 
group of 7i N k . A quasi-free state ( )q with < Q < 1 corresponds to a density matrix of the 
form 

Q ( Q Q 
t-Q \t-Q t-Q, 

The pure quasi-free states are limiting cases of this formula, their density matrices are one- 
dimensional projectors of the form 



Det(l-Q) |l©-^-©f-^-®-^-) 2 ©•••). (124) 
v ; l t-Q VI -Q t-Q) nlr J 



0©---© (p® ■■■®p\ . ©••• (125) 

fc-times 

Here P is a fc-dimensional projector acting on Hn, actually the symbol of the corresponding 
state. 

For N = 1, the algebra „4i is isomorphic to M.2- Any state p on C 2 can be mapped on a 
quasi- free state by choosing a suitable identification of A^2 with A\. Once this identification 
is fixed, the convex set of quasi-free states is equal to the set of diagonal density matrices (in 
the basis diagonalizing p). 

A2 is isomorphic to M4. One can show that any density matrix on C 4 which has the 
property that the product of its largest and smallest eigenvalues is equal to the product of the 
two others can be mapped onto a quasi-free state. Once this identification fixed, the convex 
hull of the set of quasi-free states equals 

{Ai © A2 p © A3 : (Ai, A2, A3) a probability vector and p a 2D density matrix } . (126) 
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This is exactly the set of block diagonal density matrices in a decomposition C 4 = C © C 2 © C. 

For N = 3, the convex hull of quasi-free states is also seen to be the block diagonal density 
matrices in a decomposition C 8 = C © C 3 © C 3 © C. 

For N = 4, and higher dimensions, things change. It turns out that the convex hull of 
quasi- free states is now a strict subset of the block diagonal density matrices in a decomposition 
C 16 = C © C 4 © C 6 © C 4 © C. The pure quasi-free states with support in the C 6 term correspond 
to particular ID projectors, namely restrictions of projectors the form P (g> P with P a 2D 
projector on C 4 to the antisymmetric subspace of C 4 © C 4 ~ C 6 . Such P define a 8D real 
manifold, while the ID projectors on C 6 form a 10D real manifold. Nevertheless the convex 
hull of the P © P has a non-zero volume. In fact any density matrix on C 6 may be obtained 
by a linear, generally not convex, combination of pure quasi- free states supported in C 6 . 
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